Exactly solvable model with two conductor-insulator transitions driven by impurities 
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We present an exact analysis of two conductor-insulator transitions in the random graph model. 
The average connectivity is related to the concentration of impurities. The adjacency matrix of 
a large random graph is used as a hopping Hamiltonian. Its spectrum has a delta peak at zero 
energy. Our analysis is based on an explicit expression for the height of this peak, and a detailed 
description of the localized eigenvectors and of their contribution to the peak. Starting from the 
low connectivity (high impurity density) regime, one encounters an insulator-conductor transition 
for average connectivity 1.421529 ■ ■ ■ and a conductor-insulator transition for average connectivity 
3.154985 • • •. We explain the spectral singularity at average connectivity e = 2.718281 • • ■ and relate 
it to another enumerative problem in random graph theory, the minimal vertex cover problem. 

PACS numbers: 71.30.+h, 02.50.Cw, 64.60.-i, 72.15.Rn 



Random graphs have motivated a lot of work both in 
mathematics and in physics. In the random graph model, 
N points numbered 1, 2, • ■ • , N are used as vertices. A 
pair of (distinct) vertices {i,j} is connected by an edge 
with probability p and the edges are independent. The 
adjacency matrix of the graph is the symmetric matrix 
H with matrix element H%j = 1 if vertices i and j are 
connected by an edge and zero otherwise. The average 
connectivity (i.e. average number of neighbors of a given 
vertex) is a = pN. 

An interesting asymptotic regime emerges when a is 
kept fixed as N goes to infinity. The connectivity serves 
as a parameter, and several phase transitions can be ob- 
served. 

According to the seminal papers on the subject [jjj, for 
a < 1 all connected components are finite, and only trees 
contribute to the extensive (i.e. proportional to N) quan- 
tities but for a > 1, a finite fraction of the points lies in 
a single connected component, the giant component. So 
there is a second order (classical) percolation transition 
at a = 1. 

The adjacency matrix H can be used as a Hamilto- 
nian that describes hopping of electrons from one site to 
another if the two are connected by an edge. A large 
average connectivity means small concentration of impu- 
rities and vice-versa. 

The spectrum of H is relevant to many problems in 
physics and has been investigated by many authors. It 
contains Q an infinity of delta peaks for any a and also 
a continuous component for large enough a. It has been 
argued (^] that for a ~ 1.4, a quantum percolation tran- 
sition occurs. This means that the structure of eigen- 
vectors changes : below this value, all eigenvectors are 
localized, but above, some eigenvectors occupy a finite 
fraction of the system. It is believed || that the con- 
tinuous component in the spectrum appears at the same 
threshold. Also an anomaly in the spectrum near the 
energy for a ~ 2.7 has been noticed in 

The main results presented in this letter are an analyt- 



ical expression for z(a), the average height of the delta 
peak at zero energy in the spectrum of H in the thermo- 
dynamic limit, and the identification of the corresponding 
eigenvectors, leading to precise predictions for two quan- 
tum percolation transitions in this model : a derealiza- 
tion (insulator-conductor) transition at ad = 1.421529 • • ■ 
and a relocalization (conductor-insulator) transition Q| 
at a r = 3.154985 • • •. Surprisingly, z(a) is analytic at the 
classical and quantum percolation transitions and non- 
analytic only at a = e. 

The random graph model displays some mean field 
D = oo features because the neighbors of a site are cho- 
sen with equal probability among all the other points, 
with the consequence that loops are large. This is in 
strong contrast with real materials where some underly- 
ing geometry (e.g. a lattice) restricts the possible neigh- 
bors of a point, even before impurities destroy bonds. 
But once the random graph is chosen, each point has a 
neighbors on average so it interacts effectively with only 
a finite number of points, and this is far from mean field. 
In the conducting phase, the delocahzed states live on 
the giant component, whose effective average connectiv- 
ity varies between 2.092917 • • • and 3.312453 • • •. Values 
of this order can be achieved in real two or three dimen- 
sional disordered systems, and this suggests that anal- 
ogous insulator-conductor-insulator transitions might be 
observed in real materials. Indeed, the random graph 
model gives a rich phase diagram with a single parameter, 
and some real systems should exhibit at least the same 
complexity. On the other hand, we do not expect that 
the precise values of a r and ay will appear in real mate- 
rials. Whether the critical exponents we obtain (which 
look very much like mean field exponents) remain valid 
above some upper critical dimension, and whether this 
covers some physical situation is still unclear. 

The explicit formulae below will involve five functions 
which are analytic for real positive a except maybe at 
three special values. These functions are generating func- 
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tions related to the enumeration of various types of trees. 
In particular, they all have the same small a expansion. 
We start with a short description of these functions. 

The most basic generating function for tree enumera- 
tion is the Lambert function W(a), which is analytic on 
the real positive axis and satisfies a fixed point equation, 
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(1) 



and n n is 



This implies that W = -£„>i (-a) n ^ 
the number of labeled rooted trees on n vertices. 

However, as a fixed point, W is unstable for a > e — 
exp(l), and there is a stable periodic orbit of length 2. 
This leads us to introduce the real functions A(a) and 
B{a) solving the symmetric system 



A = 



B = ae 



(2) 



and satisfying the condition that A < B for a > e. For 
< a < e, there is only one solution to the system 
(so no condition is needed in this range of a's), namely 
A(a) = B(a) — W(a). Thus e is a special value where 
a branch point occurs. We shall see later that A and 
B can also be viewed as generating functions for certain 
bicolored trees, with —a as weight for vertices of one 
color and — ae~~ a for vertices of the other color. This 
interpretation is relevant at large a. 

We shall need one more function, A* (a) which is the 
smallest solution of the equation 



A* = Ae a{e 



(3) 



Clearly A* — A is a solution, but in general not the 
smallest. The en's for which the solutions of (^) exhibit a 
branch point are such that Aa = 1 which leads to 



2 log a — ae 



— l/a 



(4) 



This equation has two solutions, ad — 1.421529- • ■, and 
a r = 3.154985 In fact, ad is the solution of the 
simpler equation 2a log a — 1. In the interval ]ay, a r [, 
A* < A, but outside this interval A* = A. The func- 
tion A* , which coincides with A for large and small a, 
and thus shares the same combinatorial interpretations in 
these regimes, appears in this discussion via the explicit 
evaluation of the sum (|l2]). 

Replacing A's by B's in (^) leads to another function 
B*(a). The functions W, A, B, A* and B* are plotted 
in Fig. 1. 

We now argue that z(a), the height of the delta peak at 
the eigenvalue zero in the spectrum of an infinite random 
adjacency matrix, is given by the explicit formula 



z(a) = -1 + 



A + B + AB 



(5) 



What we have proved rigorously B is that (g) is true 
for a < e and that the r.h.s of (||) is an upper bound for 
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FIG. 1. The special functions W , A, B, A* and B* with 
their domains of coincidence. The branchings occur at the spe- 
cial values oid,e and a r - 



z(a) for a > e. Here is the idea. Let Z(G) denote the 
number of zero eigenvalues of the adjacency matrix of 
the graph G. Suppose G has a vertex v with exactly one 
neighbor, say v' (so v is a leaf of G). Leaf removal consist 
in removing from G the vertices v and v' and all edges 
that meet those two vertices, leading to a new graph G' . 
In H it is shown that Z{G) = Z{G'). Iteration of leaf 
removal leads to a graph with consists of say I isolated 
points and a subgraph C (called the core of G) without 
leaves or isolated points. Then Z(G) = I + Z(C). In §] 
it is proved that if G is a random graph of size N and 
average connectivity a, the size of C is o(N) if Jt| and 
only if a < e and that I/N = -1 + A + B + AB + (1) for 
large N. But almost by definition, Z(G)/N = z(a)+o(l), 



and this completes the proof. 

Then z(a) is singular at a 
analytic function 

Zan{a) = -1 + 



= e. It coincides with the 



2W + W 2 



(6) 



for a < e but z(a) > z an {a) for a > e. 

We believe that (||) is true as an equality even for a > e 
for two reasons. First, the Monte Carlo simulations in 
Fig. 2 are in perfect agreement with Qq) for all a's. Sec- 
ond, from combinatorial arguments |2j,|6[, one can show 
that z(a) can be expressed formally as 



*(«) = £ 



i 



Aut(T) 



(-a) 



E{T) 



(T) 



where the sum over T is over isomorphism classes of bicol- 
ored (say brown and green) trees with at least one green 
vertex, AutiT) is the size of the automorphism group of 
T and E(T) the number of edges of T. In fact, (g) is an 
identity of analytic functions for a < 1/e. For larger a's, 
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the sum in (0) is divergent, but a natural partial resum- 
mation turns it into a convergent expression as follows. 
Any bicolored tree can be obtained from a bicolored tree 
with only green leaves by appending brown leaves at the 
green vertices, so we can explicitly sum over all brown 
leaves in the sum over bicolored trees. This yields an 
identity valid at small a: 



z{a) 



V Aut{T) 



(-a) 



VB(T) a V G (T)-l e ~V G (T) a 



(8) 



where the sum over T is over bicolored trees with only 
green leaves, and where V B (T) and Vq(T) are, respec- 
tively, the number of brown and green vertices of T. 
Regularizing the sum by restricting T to have diameter 
at most d leads to a recursive sequence depending on d 
whith limit the r.h.s of (|5|). Moreover, the sum in (||) can 
be regrouped as 



v n-l 



e~ na S n _i(-na) 



n>l 



where the Stirling polynomials S n arc defined by: 



4-n 



(9) 



(10) 



n>0 



The sum in (^|) converges to the r.h.s. of (5) for all a's, 
giving further evidence for the validity of ( 5 ) . 

The appearance of bicolored trees with only green 
leaves is no hazard. Consider the following pattern : a 
subgraph of the random graph which is a finite tree T 
and which is such that one of its two bicolorings has only 
green leaves and the green vertices share no edges with 
the complement of T in the random graph. One can show 
that the frequency of apparition of a maximal (that is, 
not contained in a larger tree with the same properties) 
such tree in the random graph is [0] 
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Aut(T) 



B V B (T) a Vfe(T)-l e - 



V G {T)a 



and that exactly Vg(T) — Vb(T) eigenvectors with eigen- 
value of the adjacency matrix of the random graph are 
localized on T (in fact, they are even localized on the 
green vertices of T) . The sum 

Zloc{a) = £ yG(: 2~^ (T) g^(r) a v g (r)-i e -v s (r)a 

(12) 

of these non-negative contributions gives a lower bound 
for z(a). Note the striking similarity between the sums 
defining z(a) and z; oc (a). Explicitly, 



zioc(a) 



Be' 



A* + B + A*B 



(13) 



and one checks that 

Ziocipi) — z(a) for < a < ad and a > a r , 
zioc(a) < z(a) for ad < a < a r . 



(14) 



This is our second important result : for < a < 
ad and a > a r , all the (extensive) contributions to the 
kernel of random adjacency matrices come from vectors 
localized on finite bicolored trees, with green leaves only, 
attached to the rest of the random graph only by brown 
vertices. However, this is not true for ad < a < a r , 
where zi oc {a) < z{a). 

Moreover, at the "critical" values ad and a r , the distri- 
bution of V B (T), V g (T), or even V G (T)-V B (T) becomes 
large : the second moment diverges. For example, 



(V G (T)-V B (T)) = 



l + 2a d - a 3 d 
a 3 , - a 2 , 



= 1.139353- •• (15) 



at connectivity a = ad but 



((V G (T)-V B (T)f) 



3a 3 d + <4 



'iad + 1 



( 2a d- a d- a d)( a d-a) 



(16) 

when a approaches ad from below. The presence of un- 
bounded fluctuations indicates that infinitely extended 
objects are responsible for the difference between zi oc (a) 
and z(a) in the range ad < a < a r . Thus we infer that at 
zero energy, the eigenvectors of the hopping Hamiltonian 
exhibit a delocalization (insulator-conductor) transition 
at ad and a relocalization (conductor-insulator) transi- 
tion at a r . The prediction of delocalization at ad is in 
numerical agreement with Monte-Carlo simulations ||. 
It turns out that the same delocalization value was al- 
ready found by Harris in a looplcss model of ran- 
dom Bethe trees. This is surprising because at ad, the 
random graph has already (many) loops. 
Note the analog of @ for zi oc (a): 



( H ) z loc (a) 



n>l 



71— 1 a —7\,0L 



n(l + B)S n - 1 {nB)-S n (nB) 



(17) 



Our main results are summarized in Fig. 2. 



As noticed before, the expressions for z[a), zi oc (a) 
and their difference show a singularity at a = e, due 
to the fact that the frequency of certain patterns is 
non-analytic. It is tempting to relate this singularity 
to the flattening of the spectral distribution near the 
zero eigenvalue for a ~ 2.7 observed in M. In fact, 
a transition in properties of random graphs at a = e 
has already been observed in a different context, mini- 
mal vertex covers ]Tc| jTTj ] . A vertex cover of a graph is 
a subset of the vertices containing at least one extrem- 
ity of every edge of the graph. An analogous concept 
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FIG. 2. Summary of the results. The functions z and zi oc 
correspond to full lines, whereas z an corresponds to the dashed 
line. The inset magnifies the region 2.5 < a < 3.5. The 
height of the delta peak z is to be compared with Monte-Carlo 
simulations, represented as black dots. The error bars are of 
order 2.10 -3 , much smaller than the symbols on the main plot, 
and of the order of the size of the symbols in the inset. The 
negative values of z an are not represented. 



is that of an edge disjoint system, i.e. a subset of the 
edges such that no two edges in the subset have a ver- 
tex in common. Given a graph G on N vertices, let 
X(G) be the minimal size of a vertex cover and Y{G) 
be the maximal size of an edge disjoint system. It is 
clear that X(G) < N and 2Y(G) < N. Those quan- 
tities share a very important property with Z(G), they 
behave simply under leaf removal || : X(G') = X(G) — 1, 
Y{G') = Y{G) - 1 whereas Z{G') = Z(G). For general 
graphs, X(G) > Y(G) and Z{G) >N- 2X{G). Exam- 
ples show that Z[G) — N + 2Y(G) can have any sign. 

But if G is a random graph with a < e, leaf removal 
leaves a core of size o(N) so if x(a), y(a) and z{a) de- 
note the limits of the averages of X(G)/N, Y(G) /N and 
Z(G)/N when N -> oo we get 
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x(a) = y(ct) = — (1 — z(a)) for a < e. 



(18) 



This is the formula for x{a) when a < e obtained 
by Hartmann and Weigt via a replica symmetric ansatz 
0,0. However, the relation ( p"8| ) has to break down at 
some a > e, because for large a, z{a) is exponentially 
small in a, whereas a rigorous result |l2|| predicts that 

2 1 

x(a) = 1 (loga - logloga - log2 + 1) + o(-). (19) 

a a 

Simulations seem to indicate that the simple relation be- 
tween x(a) and z(a) breaks down at e, together with 
replica symmetry stability. Our interpretation is that 
the remnant after leaf removal is a complicated graph 



with a size of order N which can be described only by a 
more refined replica asymmetric order parameter. 

In this paper, we have given an exact description of 
the size of the kernel of the hopping Hamiltonian asso- 
ciated to a large random graph and of the structure of 
the corresponding eigenvectors, leading to predictions of 
a derealization and relocalization transition at ad and 
a r . We have also explained the singularities affecting all 
spectral and several combinatorial properties at a = e. 
A surprising fact is that all such properties are regular 
at the (classical) percolation transition at a = 1. Pos- 
sibilities to extend these results to more realistic models 
remain to be explored. 
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